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Abstract 

It has been speculated that quantum gravity might induce a "foamy" space-time structure at 
small scales, randomly perturbing the propagation phases of free-streaming particles (such as kaons, 
neutrons, or neutrinos). Particle interferometry might then reveal non-standard decoherence ef- 
fects, in addition to standard ones (due to, e.g., finite source size and detector resolution.) In this 
work we discuss the phenomenology of such non-standard effects in the propagation of electron 
neutrinos in the Sun and in the long-baseline reactor experiment KamLAND, which jointly provide 
us with the best available probes of decoherence at neutrino energies E ~ few MeV. In the solar 
neutrino case, by means of a perturbative approach, decoherence is shown to modify the standard 
(adiabatic) propagation in matter through a calculable damping factor. By assuming a power-law 
dependence of decoherence effects in the energy domain (E n with n = 0, ±1, ±2), theoretical 
predictions for two-family neutrino mixing are compared with the data and discussed. We find 
that neither solar nor KamLAND data show evidence in favor of non-standard decoherence effects, 
whose characteristic parameter 70 can thus be significantly constrained. In the "Lorentz-invariant" 
case n = —1, we obtain the upper limit 70 < 0.78 x 10 -26 GeV at 95% C.L. In the specific case 
n = —2, the constraints can also be interpreted as bounds on possible matter density fluctuations 
in the Sun, which we improve by a factor of ~ 2 with respect to previous analyses. 

PACS numbers: 14.60.Pq, 26.65. +t, 03.65.Yz, 04.60.-m 
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I. INTRODUCTION 



Although a satisfactory theory of quantum gravity is still elusive, it has been speculated 
that it should eventually entail violations of basic quantum mechanics, including the sponta- 
neous evolution of pure states into mixed (decoherent) states [l[ through unavoidable interac- 
tions with a pervasive and "foamy" space-time fabric at the Planck scale 0] ■ The pioneering 
paper showed that such hypothetical source of decoherence might become manifest in os- 
cillating systems which propagate over macroscopical distances, through additional smearing 
effects in the observable interferometric pattern (besides the usual smearing effects due, e.g., 
to the finite source size and the detector resolution). However, lacking an "ab initio" theory 
of quantum gravity decoherence, its effects can only be parameterized in a model-dependent 
(and somewhat arbitrary) way. Searches with neutral kaon oscillations H, B|, neutron 



have 



interferometry 0, 0] and, more recently, neutrino oscillations [1, 0, 10, [ill, 112 
found no evidence for such effects so far, and have placed bounds on model parameters. 

Quantum gravity effects in neutrino systems have been investigated with increasing atten- 
tion in the last decade, as a result of the evidence for neutrino flavor oscillations. Early at- 
tempts tried to interpret the solar neutrino puzzle 0, 0, [l(| or the atmospheric v anomaly [l(| 
in terms of decoherence. After the first convincing evidence for atmospheric neutrino oscil- 
lations 15], a quantitative analysis was performed in [llj], considering possible decoherence 
effects in the — > channel (see also [12]). The phenomenology of terrestrial neutrino 
experiments has also been investigated in [l3l. 14j . More recently, prospective studies have 
focused on decoherence effects in high energy neutrinos [HI, H B, S, Hj, observable in 
next generation neutrino telescopes. Furthermore, more formal aspects of quantum gravity 
decoherence in neutrino systems have been developed [2l|, 22, [23|, 24, 25, 26, 27]. This is 
only a fraction of the related literature, which testifies the wide and increasing interest in 
the subject. 

Despite this interest, to our knowledge such decoherence effects have not been systemat- 
ically investigated in the light of the solar neutrino experiments performed in the last few 



years. Solar neutrino oscillations [28( dominated by matter effects 29l. 30| are currently well 



established by solar neutrino experiments |3ll. |32|. l33l. 134 135. |40| . |41| and have 



been independently confirmed by the long-baseline reactor experiment KamLAND [42], |43 
The striking agreement between solar and KamLAND results determines a unique solution 
in the mass-mixing parameter space [the so-called Large Mixing Angle (LMA) solution, see 



e.g. |44j, |45j], provides indirect evidence for matter effects with standard amplitude [46j, and 
generally (although not always [13]) implies that additional, non-standard physics effects 
may play only a subleading role, if any. In particular, the KamLAND collaboration has 
exploited the observation of half oscillation cycle in the energy spectrum [43J to exclude 
decoherence as a dominant explanation of their data. 

The main purpose of this paper is then to study decoherence as a subdominant effect 
in solar and KamLAND neutrino oscillations. Modifications of the standard oscillation 
formulae in the presence of decoherence, and qualitative bounds on decoherence parameters, 
are discussed in Sec. II and III for KamLAND and solar neutrinos, respectively. Quantitative 
bounds on subdominant decoherence effects from a joint analysis of solar and KamLAND 
data are studied in Sec. IV. Implications for decoherence induced by matter fluctuations in 
the Sun are discussed in Sec. V. The main results are finally summarized in Sec. VI. The 
solar neutrino flavor evolution in the presence of standard matter effects plus non-standard 
decoherence is discussed in a technical Appendix. 
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II. OSCILLATIONS WITH(OUT) DECOHERENCE IN KAMLAND 



Here and in the following, we assume the standard notation [48| for neutrino mixing, and 
set the small mixing angle #13 to zero for the sake of simplicity. For 6^3 = 0, oscillations in 
the v e —>■ v e channel probed by long-baseline reactor (KamLAND) and by solar neutrinos 
are driven by only two parameters: the mixing angle #12 and the neutrino squared mass 
difference 5m 2 = m\ — m\. In particular, the standard v e survival probability over a baseline 
L in KamLAND reads: 

P ee = l-isin 2 20 12 (l-cos^)) . (1) 

In the presence of additional decoherence effects, the oscillating factor is exponentially 
suppressed, as shown in LJJ for the atmospheric — > channel. By changing the ap- 



propriate parameters for the KamLAND v e — > v e channel, the results of lead to the 
following modification of the previous equation, 

Pee = 1 - \ sin 2 20 12 - e-* cos \ , (2) 

where the dimensional parameter 7 represents the inverse of the decoherence length after 
which the neutrino system gets mixed. 1 Equation (j2J) includes the limiting cases of pure 
oscillations (7 = and 5m 2 7^ 0) and of pure decoherence (7 7^ and Sm 2 = 0). 

Unfortunately, lacking a fundamental theory for quantum gravity, the dependence of 7 
on the underlying dynamical and kinematical parameters (most notably the neutrino energy 
E) is unknown. Following common practice, such ignorance is parameterized in a power-law 
form 

7 = 7o(#V, (3) 



vPo/ 

where E is an arbitrary pivot energy scale, which we set as E = 1 GeV in order to facilitate 
the comparison with limits on 7-parameters investigated in other contexts (as reviewed, e.g., 
in [i^]). We shall consider only five possible integer exponents, 

n = 0, ±1, ±2, (4) 

which include the following cases of interest: The "energy independent" case (n = 0); the 
"Lorentz invariant" case (n = —1); the case n = +2 that can arise in some D-brane 
or quantum-gravity models, in which 7 ~ 0(P 2 /M Planck ) ~ 10~ 19 GeV is expected (see, 



e.g. |50j); and the case where decoherence might be induced by "matter density fluctuations" 
rather than by quantum gravity {n = —2, see Sec. V). 

As previously remarked, the KamLAND collaboration [43[ (see also [HI]) has ruled out 
pure decoherence in the Lorentz invariant case (n = —1). In our statistical x 2 analysis, we 
also find that this case is rejected at 3.6a (i.e., A% 2 = 13 with respect to pure oscillations). 
In addition, we find that the other exponents in Eq. (j4j) are also rejected at > 3cr for the pure 
decoherence case. Therefore, decoherence effects can only be subdominant in KamLAND, 
namely 

< 1 . (5) 



1 Units: [7] = 1/length = energy. Conversion factor: (1 km) 1 = 1.97 x 10 19 GeV. 
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For typical KamLAND neutrino energies (E ~ few MeV) and baselines (L ~ 2x 10 2 km), the 
above inequality implies upper bounds on 70, which range from 70 <C 1CT 26 GeV {n = —2) 
to 70 <C 10~ 16 GeV {n = +2). We do not refine the analysis of such bounds (placed by 
KamLAND alone), since they are superseded by solar data constraints, as shown in the next 
Section. 

III. OSCILLATIONS WITH(OUT) DECOHERENCE EFFECTS IN SOLAR NEU- 
TRINOS 

The survival probability describing standard adiabatic v e transitions in the solar matter 
is given by the simple formula (up to small Earth matter effects) 

P® = \ (l + cos2Mr ) cos2# 12 ) , (6) 

where ^12(^0) is the energy- dependent effective mixing angle in matter at the production 
radius r (see, e.g., [52j and references therein). 

In the presence of non-standard decoherence effects, we find that the energy dependent 
term is modulated by an exponential factor, 

if e = \ (l + e"^ e cos20 12 (r o ) cos2# 12 ) , (7) 
where R Q = 6.96 x 10 5 km is the Sun radius, while 7 is defined as 

7© = 7o9n(E) , (8) 

where the dimensionless function g n {E) embeds, besides the power-law dependence E n , also 
the information about the solar density profile [which is instead absent in Eq. (6)]. The 
reader is referred to the Appendix for a derivation of Eq. and for details about the 
function g n (E). 

Equation ([7]) includes the subcase of pure oscillations (70 = and 5m 2 7^ 0), but not the 
subcase of pure decoherence (70 7^ and 5m 2 = 0), since the limit 5m 2 — > would entail 
strongly nonadiabatic transitions, and thus a breakdown of the adiabatic approximation 
assumed above. However, as noted in the previous section, the KamLAND data exclude the 
limit 5m 2 — > 0; moreover, they require 5m 2 values which are high enough to guarantee the 
validity of the adiabatic approximations, with or without subleading decoherence effects (as 
we have numerically verified). Therefore, for solar neutrinos, the only phenomenologically 
relevant cases are those including oscillations plus decoherence. 2 

We have analyzed all the available solar neutrino data with (5m 2 , #12, 70) taken as free 
parameters. It turns out that, despite the allowance for an extra degree of freedom (70), the 
data always prefer the pure oscillations (70 = 0) as best fit, independently of the power-law 
index in Eq. (j4j). Since there are no indications in favor of decoherence effects, the exponent 
in Eq. ([7]) is expected to be small, 

70^0 < 1 ■ (9) 



2 For the sake of curiosity, we have anyway calculated P® for the pure decoherence case, by numerically 
solving the neutrino evolution equations (discussed in the Appendix) for 8m 2 — and 7 7^ 0. We always 
find P® > 1/2, which is forbidden by 8 B solar neutrino data [4q . 
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TABLE I: Upper limits on the decoherence parameter 70 obtained for different values of n from a 
global fit to solar and KamLAND data, after marginalization of the mass-mixing parameters. The 
limits refer to 95% C.L. (i.e., 2a, or A% 2 = 4). 



n 70 (GeV) 

-2 < 0.81 x 10- 

-1 < 0.78 x 10~ 

< 0.67 x 10~ : 

+1 < 0.58 x 10- 

+2 < 0.47 x 10~ 



For typical solar neutrino energies E ~ 10 MeV, it turns out that g n {E) ~ 0.2 x 10 -2 " (see 
the Appendix), and the above inequality can be translated into upper bounds on 70, which 
range from 70 < 10~ 28 GeV (for n = —2) to 70 10~ 20 GeV (for n = +2). Such bounds 
are two to four orders of magnitude stronger than those placed by KamLAND alone (see 
the end of the previous Section). A useful complementarity then emerges between solar and 
KamLAND data in joint fits: The former dominate the constraints on decoherence effects, 
while the latter fix the mass-mixing parameters independently of (negligible) decoherence 
effects. 



IV. COMBINATION OF SOLAR AND KAMLAND DATA: RESULTS AND DIS- 
CUSSION 

We have performed a joint analysis of solar and KamLAND data 3 in the (5m 2 , sin 2 # 12 , 70) 
parameter space for the five power-law exponents n = 0, ±1, ±2. The main results are: (i) 
70 = is always preferred at best fit, i.e., there is no indication in favor of decoherence effects; 
(ii) the best fit values and the marginalized bounds for (5m 2 , sin 2 #12) do not appreciably 
change from those obtained in the pure oscillation case, namely, 5m 2 = (7.92 ± 0.71) x 10 -5 
eV 2 and sin 2 6>i 2 = (O.SM+q^) at ±2a [44]; (hi) significant upper bounds can be set on 
the decoherence parameter 70 . Our limits on 70 are given numerically in Table I (at the 2a 
level, Ax 2 = 4) and graphically in Fig. 1 (at 2a and 3a level). Such limits are consistent 
with those discussed qualitatively after Eq. §§§ in the previous Section. 

Figure 1 clearly shows that the bounds on 70 scale with n almost exactly as a power law, 
changing by about two decades for |An| = 1. The reason is that the bounds are dominated 
by solar neutrino data, and in particular by data probing the 8 B neutrinos in a relatively 
narrow energy range around E ~ 10 MeV; the power-law dependence assumed in Eq. ([3]) 
and embedded in the function g n (E) then implies that the parameter 70 scales roughly as 
(E/E )~ n ~ 10 2n . 

Although the case of no decoherence (70 = 0) is preferred, it makes sense to ask what 
one should observe for decoherence effects as large as currently allowed by the data at, say, 
the 2a level. Figure 2 compares the P® energy profile for the cases of pure oscillations (left 
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The details of the data set, of the solar model used [53( and of the statistical x 2 analysis have been 



reported in 44j and are not repeated here. 
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panel) and of oscillations plus decoherence (right panel), where 70 is taken equal to the upper 
bound at 2a, as taken from Table I. 4 It can be seen that decoherence effects, to some extent, 
mimic the effects of a larger mixing. For instance, the curve with n = in the right panel is 
not much different from the curve at sin 2 # 12 = 0.371 (upper 2a value) in the left panel. As 
a consequence, one expects some degeneracy between the parameters 70 and sin 2 6*12 when 
fitting the data. The degeneracy is only partial however, because decoherence effects can 
significantly change both the shape and the slope of the energy profile within the current 
2a bounds, as evident in the right panel. Therefore, future measurements of the (currently 
not well constrained) solar neutrino energy spectrum will provide further important probes 
of decoherence effects. 

The variation of P® due to subdominant decoherence effects [see Eqs. ([6]) and ([7])] is 
given, in first approximation, by 

AP® ~ -~ 7o P cos 20i 2 cos 29 12 , (10) 

and changes sign with cos2#i2. As the energy increases, the value of cos2#i2 changes from 
cos 812 > (low-energy, vacuum-dominated regime) to —1 (high-energy, matter-dominated 
regime), the transition being located around 2 MeV for the 8 B neutrino curves shown in 
Fig. 2. This fact explains the general increase of P® for E > 2 MeV in the right panel of 
Fig. 2. More detailed features depend instead on the energy behavior of the function g n (E), 
which modulates decoherence effects (see the Appendix). In general, g n (E) grows rapidly 
with increasing energy for n > (which explains the high-energy upturn of the curves with 
n = +1 and n = +2), while it vanishes with decreasing energy for all n 7^ —2 (which explains 
the low-energy equality of all curves but for n = —2). 5 The "bunching" of the curves around 
E ~ 10 MeV in the right panel of Fig. 2 is in part a data selection effect, since this energy 
region is strongly constrained by precise 8 B neutrino data. Further spectral 8 B data will 
be very useful to constrain the slope of the energy spectrum and thus also the sign of the 
power-law index n. 

Figure 3 illustrates the partial degeneracy between decoherence effects and mixing angle, 
as a shift in the allowed regions for fixed 70 7^ in three representative cases (from left 
to right, n = —2, 0, +2). In each panel, the thin dotted curves enclose the mass-mixing 
parameter regions allowed at 2a by the standard oscillation fit of solar data (larger region) 
and by solar plus KamLAND data (smaller region). The thick solid curves refer to the 
same data, but fixing a priori the decoherence parameter 70 at the 2a upper limit value 
in Table I. In all cases in Fig. 3, the curves with 70 7^ are shifted to lower values of the 
mixing angle, as compared to pure oscillations; this means that decoherence effects can be 
partly traded for a smaller value of the mixing angle in solar neutrino oscillations. Therefore, 
should future solar neutrino data prefer smaller (larger) values of sin 2 #12 with respect to 
KamLAND data, there would be more (less) room for possible subdominant decoherence 
effects. As already remarked, the degeneracy between 70 and sin 2 8\ 2 is only partial, and 
future neutrino spectroscopy will provide a further handle to break it, should decoherence 
effects (if any) be found. 



4 For definiteness, Fig. 2 shows the daytime probability of 8 B neutrinos, averaged over their production 
region in the Sun. 

5 Constraints on the specific case n — — 2 might thus benefit of sub-MeV solar neutrino observations in 
Borexino 54 1. 
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We conclude this section by confronting the bounds in Table I with those derivable from 
the analysis of atmospheric neutrino data (which, by themselves, exclude pure decoherence, 
at least in the n = —1 case [551] ). In principle, a direct comparison is not possible, since the 
7o parameter introduced here for the solar v e — > v e channel does not need to be the same as 
for the atmospheric — > channel. However, if the 70's for these two channels are assumed 
to be roughly equal in size, then it is easy to realize that solar+KamLAND neutrinos set 
stronger (weaker) bounds than atmospheric neutrinos for n < (n > 0), as a consequence 
of the different neutrino energy range probed. In fact, due to the assumed power-law energy 
dependence of decoherence effects, negative (positive) values of n are best probed by low- 
energy solar (high-energy atmospheric) neutrino experiments. For the intermediate case 
(n = 0), it turns out that matter effects render solar neutrinos more sensitive to 70 than 
atmospheric neutrinos. Just to make specific numerical examples: for n = (—1, 0, +2) one 
roughly gets the bounds 70 ^ (0.7 x 10~ 21 , 0.4 x 10~ 22 , 0.9 x 10~ 27 ) GeV from atmospheric 
[Til (plus accelerator 12]) neutrino data, to be compared with the corresponding limits from 
solar+KamLAND data from Table I, 70 < (0.78 x 10" 26 , 0.67 x 10~ 24 , 0.47 x 10" 20 ) GeV. 
Solar neutrinos clearly win over atmospheric neutrinos for n < 0. This comparison must be 
taken with a grain of salt, since it can radically change by assuming either independent 70's 
in different oscillation channels, or functional forms of j(E) different from power laws. 

Finally we observe that in the n = +2 case, motivated by some "quantum-gravity" or 



"string- inspired" models [50|, the solar+KamLAND limit on 70 is one order of magnitude 
lower than the theoretical expectation [70 ~ O(i^o/Mp lanck ) ~ 10~ 19 GeV]. In the case of 
atmospheric neutrinos, this bound is even stronger (70 < 0.9 x 10~ 27 GeV). Consequently, 
these models appear strongly disfavored, at least in the neutrino sector. 



V. RECOVERING THE CASE OF DENSITY FLUCTUATIONS IN THE SUN 



Decoherence effects in solar neutrino oscillations can be induced not only by quantum 
gravity, but also by more "prosaic" sources, such as matter density fluctuations — possibly 
induced by turbulence in the innermost re gion s of the Sun. This topic has been widely 



investigated in the literature [56l . |5j 



59l . 60, 61, 62, 111, 64, 65|, and quantitative upper 



limits have already been set [6l|, [62|, 63|, 64j by combining solar data and first KamLAND 
results. 

It turns out that stochastic density fluctuations lead (with appropriate redefinition of 
parameters) to effects which have the same functional form as those induced by quantum 
gravity in the n = —2 case. More precisely, let us consider fluctuations of the solar electron 
density N e around the average value (N e ) predicted by the standard solar model, 



N e (r) = (l + 0F(r))(N e (r)), 



(11) 



where F(r) is a random variable describing fluctuations at a given radius r, and /3 represents 
their fractional amplitude around the average. It is customary to assume a delta-correlated 
(white) noise, 

(F(r 1 )F(r 2 )) = 2r5(r 1 -r 2 ), (12) 

where r is the correlation length of fluctuations along the (~ radial) neutrino direction. 

As shown in [59| through a perturbative method (which inspired our approach to deco- 
herence in the Appendix), the effect of delta-correlated noise on adiabatic neutrino flavor 
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transitions can be embedded through an exponential damping factor as in Eq. (JTj). The func- 
tional form turns out to be the same as for the n = —2 case, provided that one makes — in 
our notation — the replacement 

fr(^X . (13) 



7o 



2E 



By using the bound 70 < 0.81 x 10 28 GeV (Table I, case n = —2) and the best-fit value 
5m 2 = 7.92 x 10 -5 eV 2 , one gets the following upper limit 



f3 2 r < 1.02 x 10~ 2 km (2a) , 



(14) 

on the parameter combination (3 2 t which is relevant 0, 63[ for density fluctuation effects 
on neutrino propagation. 



As stressed in 63 



care must be taken in extracting an upper limit on the fractional am- 
plitude (3 for fixed correlation length r. Indeed, the delta-correlated noise is an acceptable 
approximation only if the correlation length r is much smaller than the oscillation wave- 
length in matter — a condition that becomes critical for low neutrino energies. In particular, 



assuming a reference value r = 10 km as in |63j, such condition is violated for energies 
E ^$ 1 MeV. Following 63], we have thus excluded low-energy, radiochemical solar neutrino 
data form the solar+KamLAND data fit, and obtained a slightly weaker (but more reliable) 
upper bound from the analysis of 8 B neutrino data only, 



(3 2 r < 1.16 x 10~ 2 km (2a) , 



(15) 



which, for r = 10 km, translates into an upper limit on the fractional fluctuation amplitude, 

P < 3.4% (2a) . (16) 



This limit improves the previous one derived in [63] ((3 < 6.3% at 2a) by a factor of ~ 2, 
essentially as a result of the inclusion of the most recent solar and KamLAND neutrino data 
appeared after [631] . 6 



VI. SUMMARY AND CONCLUSIONS 

In this paper, we have investigated hypothetical decoherence effects (e.g., induced by 
quantum gravity) in the v e — > v e oscillation channel explored by the solar and KamLAND 
experiments. In both kinds of of experiments, decoherence effects can be embedded through 
exponential damping factors, proportional to a common parameter 70, which modulate the 
energy-dependent part of the v e survival probability. By assuming that the (unknown) 
functional form of decoherence effects is a power- law in energy {E n ), we have studied the 
phenomenological constraints on the main decoherence parameter (70) for n = 0, ±1, ±2. 
It turns out that both solar and KamLAND data do not provide indications in favor of 
decoherence effects and prefer the standard oscillation case (70 = 0) for any index n. By 
combining the two data sets, n-dependent upper bounds (dominated by solar neutrino data) 



6 We have verified that, by adopting the same (older) data set and standard solar model as used in 63j, we 
recover the same 2a upper limit, (3 < 6%. 



S 



have been derived on 70, as reported in Table I and shown in Fig. 1. In the "Lorentz- 
invariant" case n = — 1, we obtain the upper limit 70 < 0.78 x 1CT 26 GeV at 95% C.L. For 
n = —2, the results can also be interpreted as limits on the amplitude of possible (delta- 
correlated) density fluctuations in the Sun, which we improve by a factor of two [Eqs. (1T5]) 
and (fTB"j) ] with respect to previous bounds. 

Further progress might come from a better determination of the energy profile of solar 
neutrino flavor transitions as well as from more precise measurements of sin 2 # 12 (which 
is partly degenerate with 70), attainable with KamLAND and future long-baseline reactor 



neutrino experiments 66 



APPENDIX: DECOHERENCE AND MATTER EFFECTS IN SOLAR NEUTRI- 
NOS 



In this section we discuss a perturbative calculation of decoherence effects for solar neutri- 
nos, where matter effects are known to be relevant. The approach, inspired by the work 59J, 
draws on the formalism and the notation introduced in for the case of decoherence in 
the — > channel, here adapted to the v e — > v e channel. In the following, the the notation 
is made more compact by setting 9 = 9 12 , c 2 e = cos2# 12 , s 2 g = sin 2#i 2 , P ee = -P e ® etc. 

Decoherence effects in the flavor evolution of the (z/ e , v a ) system (where a = /x, r) along 
the space coordinate r(~ t) 7 can be described in terms of the neutrino density matrix, 
obeying a modified master equation of the form 67 



dp 
dr 



JL = -i[H v + H n (T),p\- 1 [D,[D,p\\ 



(17) 



where H v (H m ) is the "vacuum" (matter) Hamiltonian, and the operator D embeds de- 
coherence effects with amplitude 7, parameterized as in Eq. 7 = •y Q (E/E ) n with 
Eq = 1 GeV. While unitarity is preserved (i.e., Trp(r) = 1), coherence is lost in the propa- 
gation (J^Trp 2 < 0). Equation (fTTj) satisfies the conditions of complete positivity [68j and 



non-decreasing entropy in the v system evolution [6J 
In the flavor basis the standard oscillation terms read 



H v 



~U e a 3 U t 
V(r) 

7T~ °"3 



-c 2 e s 2d 
S20 c 29 



(18) 
(19) 



where 03 is the third Pauli matrix, k = 5m 2 /2E is the vacuum wavenumber, and V(r) = 
\^2GFN e (r) is the interaction potential in matter. 

As in [ill ], we assume energy conservation for evolution in vacuum, i.e., Ti[H v p(r)] = con- 
stant. This condition is satisfied if [H v , D] — 0, 70l |. namely, in a two-dimensional system, 
if D oc H v . We can thus take 



Da 



-c 2S s 29 
S29 c 2 e 



(20) 



7 Note that r does not necessarily coincide with the radial coordinate, due to the extended neutrino pro- 
duction region (which is taken into account in our analysis.) 
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without loss of generality, since any overall factor can be absorbed in 70. We also make the 
plausible assumption that Eq. (1201) is not altered for evolution in matter, since decoherence 
induced by quantum gravity is unrelated to electroweak matter effects. 

As usual, Eq. (IT7|) can be written in terms of a "polarization" vector P with components 
Pi = |Tr[p<Tj] (Bloch equation): 



dP 

dr 



where n = [s 20 , 0, -c 2e ] T , P_ 



H(r) 



[kn + V(r)e 3 ) x P - 7 Pj 
H(r)P - 7 £> e P 

(P ■ n)n, and: 





V(r) - kc 2 e 




-V{r) + kc 2e 

ks 26 



(21) 





-ks 2e 




(22) 



-2B 





C 2 gS 2e 

1 



c 2 es 2 e s 



26 



(23) 



Note that for V = (vacuum propagation), the solution of the above Bloch equation leads 
to the survival probability in Eq. d2J), (see [ll| for details.) 

The matrix 7i has eigenvalues Aq = and A± = ±k, corresponding to the eigenvectors: 



u 



'261 



C 2§ 



1 
V2 



C 2B 

±i 

S 29 



(24) 



In the above equations, a "tilde" marks effective parameters in matter: k is the oscillation 
wavenumber in matter, defined through k/k = [1 — 2Vc 2 e/k + (V/k) 2 ] 1 ^ 2 , while 9 is the 
mixing angle in matter, defined through s 2 q = ks 2 g/k and c 2 ^ = (kc 2 g — V)/k. The matrix 7i 
is diagonalized through the matrix 11(9) = [u ,u + ,u_]: TV (9) ■ H ■ TZ(9) — diag[0, +ik, — ik) . 

In the absence of decoherence effects, the adiabatic solution of Eq. ([2T]) appropriate for 
current solar neutrino phenomenology is 



P(R Q ) = K{0) ■ diag 



1, e 



+i f R@ drk(r) -i \ Rq drk(r) 
Jr . e r 



■TV(9aMro) 



(25) 



where P(ro) = T [0,0, 1] for an initial v e state, 9q is the mixing angle in matter at the 
production point r , and 9 = 9 is taken at r = R®. After averaging on the fast oscillating 
terms (a "standard" decoherence effect), one recovers the usual adiabatic formula [Eq. 
for the survival probability P ee 



Tr\p\v e )(v e \] 



1 + P 3 (R & ) 1 + c 2e c. 



29 



(26) 



Let us now treat the term —^DgP in Eq. (j2"T|) as a perturbation [59( . The corrections to 
the eigenvectors lead to variations of P ee of 0(j/k) ^ 10~ 3 (for the range of j/k allowed a 
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posteriori by the fit to solar neutrino data) and can be neglected. The corrections to the 
two eigenvalues A± can also be neglected, since they would only lead to a further damping 
of the fast oscillating terms, which are already averaged out. 8 

The first-order correction to the eigenvalue Ao (whose unperturbed value is zero) is the 
only relevant one, 

5\ = --fulV e u = -7 sin 2 2(0 - 9) . (27) 
and leads to the following correction to Eq. fl26|) : 



where 



7 



1 




e 


2 




V(r)s 2§ (r) 


Jro 


k 



dr 



(28) 



(29) 



Equation (JTj) is then recovered by setting 7 
function g n (E) as 



7o<7 n (i?) and by defining the dimensionless 





V{r)s 2§ (r) 


\rjQj Jr 


k 



1 2 



dr 
R Q 



(30) 



The function g n {E) depends mostly on the neutrino energy E and, to some extent, on 
the parameters r , 5m 2 , and sin 2 # 12 . Figure 4 shows this function as calculated for r = 0, 
Sm 2 = 7.92 x 10" 5 eV 2 , and sin 2 12 = 0.314. For E -> the function g n (E) (and the 
associated decoherence effect) vanishes, except for the case n = —2, where the factors E~ 2 
and k~ 2 cancel out and provide a finite limit (7-2(0) 7^ 0. 

We have tested the analytical Eq. (!28l) against the results of a numerical integration of 
the Bloch equation, for many representative points in the parameter space relevant for solar 
v phenomenology, and we find very good agreement (SP ee < 10 -4 ) for all values of n 7^ —2. 
Only in the case n = —2, the comparison of analytical and numerical results is slightly 
worse (but still very good, 5P ee < 10 -3 ) at the lowest detectable energies (~ 0.1 MeV), 
due to the breakdown of perturbation theory for E — > 0. For practical purposes, however, 
the modified adiabatic Eq. (128]) accurately replaces the results of numerical solutions of the 
Bloch equation for solar neutrinos. 
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Solar + KamLAND constraints on y 
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FIG. 1: Upper bounds on the decoherence parameter 70 as a function of the power-law index n, 
as obtained from a combined analysis of solar and KamLAND data. The solid and dotted curves 
refer to 2a and 3a confidence level, respectively. 
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Standard oscillations (y = 0) Oscillations + decoherence 
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FIG. 2: Energy profile of the (daytime) survival probability of 8 B neutrinos, averaged over their 
production region in the sun: Comparison of the effects produced by variations of sin 2 6 12 for 
70 = (left panel) and by 70 7^ at fixed sin 2 612 (right panel). In the left panel, sin 2 9i% is varied 
within its ±2<t limits. In the right panel, for each index n = [—2, —1,0, +1, +2] the value of 70 is 
taken equal to the corresponding 2a upper limit (reported in Table I), which in units of 10 -24 GeV 
corresponds respectively to: 0.81 x 10~ 4 (n = -2), 0.78 x 10" 2 (n = -1), 0.67 (n = 0), 0.58 x 10 2 
(n = +1), 0.47 x 10 4 (n = +2). The curve corresponding to standard oscillations (70 = 0) is also 
shown in the right panel as a guide to the eye. In all cases, 5m 2 is fixed at its best-fit value. 
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FIG. 3: Constraints on the mass-mixing parameters for 70 = (thin dotted curves) and for 70 
fixed at its 2a upper limit in Table I (thick solid curves). The three panels refer, from left to 
right, to the three cases n = —2, 0, and +2. The smaller (larger) allowed regions refer to the 
solar+KamLAND (solar only) data analysis. 
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FIG. 4: Energy profile of the auxiliary function g n (E), which modulates the exponent of the 
damping factor induced by decoherence in solar neutrino oscillations. The function is multiplied 
by 10 3n for a better graphical view. The shown function refers to a neutrino produced at the Sun 
center and to best-fit oscillation parameters; altering this choice would induce minor variations. 
See the text for details. 
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